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^ ' Abstract. To a system of second order ordinary differential equations (SODE) one can assign 

Q ' a canonical nonlinear connection that describes the geometry of the system. In this work we 

, develop a geometric setting that allows us to assign a canonical nonlinear connection also to a 

1 system of higher order ordinary differential equations (HODE). For this nonlinear connection 

' we develop its geometry, and explicitly compute all curvature components of the correspond- 

^N) ' ing Jacobi endomorphism. Using these curvature components we derive a Jacobi equation that 

describes the behavior of nearby geodesies to a HODE. We motivate the applicability of this 
. nonlinear connection using examples from the equivalence problem, the inverse problem of the 

' calculus of variations, and biharmonicity. For example, using components of the Jacobi endo- 

' morphism we express two Wuenschmann-type invariants that appear in the study of scalar third 

or fourth order ordinary differential equations. 

a 

P.- 

1. Introduction 

It is well known that to a system of second order ordinary differential equations (SODE) one 
^ ' can associate a nonlinear connection in a canonical way. This nonlinear connection further induce 

0^ , various geometric objects useful for the study of the SODE such as the dynamical covariant deriv- 

' ative, parallel transport, torsions, curvatures, and the Jacobi equation [71 HZl [551 [271 [3S] . There 

are various attempts to define a similar canonical nonlinear connection also to systems of higher 
order ordinary differential equations (HODE), [5l[6l[l8l[20l[3Tl[33l[34|. A difficulty with this task 
is that to a HODE one can associate a number of nonlinear connections, and each of these offer 
different information about the HODE, [H H H [H [H [35] . 

The purpose of this work is to develop a unifying geometric setting for studying systems of second 
and higher order ordinary differential equations. Our approach is focused on two ideas of Kosambi 
[271 128|. First, in [27], Kosambi states that a "system of differential equations can be dealt with 
' geometrically by means of the tensorial operator of differentiation, called the bi-derivative". In this 

work we propose a global expression, called the dynamical covariant derivative, for the bi-derivative 
operator introduced by Kosambi. The dynamical covariant derivative we propose, is associated to 
a system of differential equations, of second or higher order, and a nonlinear connection that is not 
fixed yet. We determine the nonlinear connection by requiring the compatibility of the dynamical 
covariant derivative with some geometric structure, introduced by the equations of variation of the 
system of differential equations. This follows the statement of Kosambi, [28], that for a system of 
HODE, "the connection cannot be determined without the recourse to the equations of variation" . 

In Section [5] we review the geometrical setting for systems of second order ordinary differential 
equations. We describe the canonical nonlinear connection [T71 124| . the dynamical covariant de- 
rivative [T0l[T4], symmetries, newtonoid vector fields [32l I40|. the Jacobi endomorphism and the 
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Jacobi equation |191 US]. The purpose of this section is to present these objects such that their 
generahzation to systems of higher order ordinary differential equations will be natural. 

In Section 13] we start with a system of HODE, identified with a vector field S, which is called 
a semispray, on the k-th order tangent bundle T^M of the configuration manifold M. In this 
part of the work we show what are the advantages of using a specific nonlinear connection if we 
want to obtain information about the first order variation and the symmetries of a system of 
higher order differential equations. This connection has been proposed in a coordinate form by 
Miron and Atanasiu in and has its origins in the work of Kosambi [2S]. In Proposition 13.81 
and Theorem 13.101 we fix this nonlinear connection, among the other ones that can be associated 
to a semispray of order k, using the following structures: the dynamical covariant derivative V, 
the Jacobi endomorphism and an algebraic structure of the set of newtonoid vector fields 
of order k. In Theorem 13.111 we use the chosen nonlinear connection to characterize the set of 
symmetries of a semispray of order k using higher order dynamical covariant derivatives and the 
Jacobi endomorphism. We provide explicit formulae for all curvature components of the Jacobi 
endomorphism. In the last part of this section we discuss the first order geodesic variation and 
obtain the corresponding Jacobi equations for a semispray of order k. 

In Section |4] we discuss the applicability of the theory developed in this paper in various fields, 
using examples from the equivalence problem, the inverse problem of the calculus of variation, and 
biharmonicity. In Subsection 14.11 we consider a second order Lagrangian L2 that represents the 
second order prolongation of a Riemannian structure (5] . Its variational semispray is a semispray 
of order 3 and the corresponding Jacobi endomorphism encodes information about the curvature 
of the Riemannian structure, and their first and second order covariant derivatives. The second 
order Lagrangian L2 has been used also in |13) to prove that biharmonic curves are variational. 
We further show that some invariants that have been associated to third or fourth order ordinary 
differential equations can be expressed in a geometrical way in terms of the curvature components of 
the Jacobi endomorphism and their dynamical covariant derivatives. In Subsection 14.21 we provide 
a geometric reformulation for one of these invariants, called the Wuenschmann invariant. The 
Wuenschmann invariant represents the obstruction for the existence of an associated conformal 
Lorentzian structure on the 3-dimensional solution space of a third order differential equation, 
PPI [53]. Also, it has been shown that the Wuenschmann invariant represents an obstruction 
for an equivalence problem of third order ordinary differential equations, |38j . In Subsection 
14.31 we provide a geometric reformulation for an invariant that was proposed by Fels in [22] for 
characterizing variational fourth order ordinary differential equations. This invariant has been also 
considered by Dridi and Neut in |21) for studying the equivalence problem of fourth order ordinary 
differential equations. 

2. Symmetries and first order variation for systems of second order differential 

equations 

In this section we reformulate some results about the geometry of a system of second order 
ordinary differential equations such that each of them has a generalization to system of higher 
order ordinary differential equations. 

2.1. Systems of SODE and geometric structures on TAI. Consider M a real, n-dimensional 
manifold, which is a topological Hausdorff space with countable base that is locally homeomorphic 
to M", and with C°°-smooth transition maps. We assume that all objects are smooth where 
defined. For a manifold M, its tangent bundle (TM, tt, M) is denoted by TM, the ring of smooth 
functions on M is denoted by C°°(M), and the set of vector fields on M is denoted by j£(M). The 
canonical submersion of the tangent bundle, tt : TM — )■ M, induces a natural foliation on TM, 
whose leaves are tangent spaces 7r~^(p) = TpM , iov p G M. Local coordinates on the base manifold 
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M are denoted by (x*), while induced coordinates on TM will be denoted by Consider 
also TM \ {0} = {{x, y) 6 TM, y ^ 0}, the tangent bundle with the zero section removed. 

Throughout the paper we will consider / a nonempty open interval in M. For a curve c : / — >■ M , 
c{t) = {x'^{t)), denote by c' : / — J- TM, c' {t) = dx'^/dt), its tangent lift. We say that curve c 

is a regular curve if c'(t) e TA/ \ {0} for all t G /. 

The tangent structure (or vertical endomorphism) is the (1, l)-type tensor field on TM defined 
as J = d/dy^ (i) dx^. The vertical distribution is the n-dimensional distribution given by V : u € 
TM VuTM = Kcr J„ = ImJu = Ker DuTt C TuTM, where DuTt denotes the tangent map of tt 
at M G TAI. This distribution is integrable being tangent to the natural foliation of the tangent 
bundle. 

In this section we pay attention to the geometry of a system of second order ordinary differential 
equations. Such a system of SODE on M can be represented using a semispray, which is a special 
vector field on TM\{0}. 

A semispray on M is a vector field S G X{TM \ {0}) such that any integral curve of S, 
7 : / — > TM \ {0}, is of the form 7 — (tt o 7)'. The reason that a semispray S is only defined 
on TM \ {0} and not on the whole tangent bundle TM is to include also semisprays induced by 
Finsler metrics [4], I43j. In induced coordinates (a;, y) for TM \ {0}, a semispray S is given by 

(2.1) S^fl-^2G^ix,y)^, 

for some functions defined on domains of {x, y). For an integral curve 7 : / — > TM \ {0} of 5, 
we say that curve c : I M , c = tt o 7 is a geodesic of S. Therefore, a regular curve c : / — A/ is 
a geodesic of S if and only if S* o c' = c". Locally, a regular curve c : I ^ M, c{t) = (.T*(t)), is a 
geodesic of 5 if and only if it satisfies the system of second order ordinary differential equations 

Thus a semispray describes systems of SODE with regular curves on M as solutions. 

A geometric structure useful for the geometry of a system of SODE is that of nonlinear con- 
nection. An n-dimensional distribution H : u £ TM \ {0} 1-^ HuTM C TuTM that is supple- 
mentary to the vertical distribution is called a nonlinear connection, or horizontal distribution. 
Therefore, for a nonlinear connection, we have the following decomposition: TuTM ~ Hu © Vu, 
for all u G TM \ {0}. We will denote by h and v the horizontal and vertical projectors that 
correspond to the above decomposition, and since h + v = Id, either of these two projectors 
characterizes the nonlinear connection H. We will denote by X^{TM \ {0}) = h {X{TM \ {0})) 
and XyiTM \ {0}) = v {X{TM \ {0})) the C°°{TM \ {0})-modules of horizontal and respectively 
vertical vector fields. 

In induced coordinates {x\ y') for TM, the vertical distribution is spanned by d/dy\ Therefore, 
the horizontal distribution is locally spanned by vector fields of the form 

Sx^ ~ \dx^ ) ~ dx^ ' 

for some functions Nj defined on domains of induced coordinates on rA/\{0}. The corresponding 
horizontal and vertical projectors can be written as 

S , , 9 - ,■ 

h=—(»dx\ v = —®5y\ 
ox ay 

The 1-forms 5y^ — dy' + Njdx^ are annihilators for the horizontal distribution. A change of induced 
local coordinates {x^ ,y'^) — {x^ {x) ,y^ {x , y)) on TM induces the following transformation rule for 
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the basis 5/5x'', d/dy^ and its dual basis dx\Sy^ 



6 9a;' 6 ^ ^ dx^ d _ dx^ ^ ^ — § i 

5x^ dxi Sx"^ ' dy^ dxi dy^ ' dx^ ' dx^ 

Therefore the components of a tensor field on TAf , with respect to these bases, will transform as 
the components of a tensor field on the base manifold M . 

Although a semispray S always induces a canonical nonlinear connection, we will not yet assume 
any relation between S and (/i, v) at this point. For a semispray S and an arbitrary nonlinear 
connection {h,v), consider the vertically valued (1, l)-type tensor field on TAI \ {0}, 

<I> = —V o Csv = V o Csh ~ V o Cs ° h, 

which will be called the Jacobi endomorphism. Here, for a (f , l)-type tensor field A, we denote 

CsA^ CsoA- AoCs, 

the Frolicher-Nijenhuis bracket of 5* and A. Locally, the Jacobi endomorphism can be expressed 
as follows 

(2.3) $ = R^^2j--S{Nj) + NlN^. 

The components i?* of the Jacobi endomorphism have been considered be Kosambi [27j as the 
second geometric invariant of the semispray. From this invariant one can determine two other 
geometric invariants, namely the curvature of the nonlinear connection and one component of the 
curvature of the Berwald connection [5]. In the Riemannian context, the components i?* of the 
Jacobi endomorphism are related to the curvature i^j-j.;, |37j . of the Levi-Civita connection as 
follows R^j{x,y) = R\^i{x)y^yK A similar formula holds true also in the Finslerian context, [i^ 
(8.15)], where i?* is called the Riemann curvature. For the Jacobi endomorphism $ in formula 
(|2.3p . we will refer to its components as to the curvature components. 

For a pair (5, (/i, w)), consider the map V : X(rM \ {0}) -J> X(rM \ {0}), given by 

(2.4) y — ho Cs o h + V o Cs o v ^ Cs + ho Csh + v o Csv 

that will be called the dynamical covariant derivative. By setting V/ ~ S{f), for / G C°°{TM \ 
{0}), using the Leibniz rule, and the requirement that V commutes with tensor contraction, we 
extend the action of V to arbitrary tensor fields and forms on TM \ {0}, see |101 Section 3.2]. For 
example, if w is a 1-form on TM \ {0}, then its dynamical covariant derivative is given by 

(2.5) {\7u){X) = S{u{X))-u{\7X). 

For a (1, l)-type tensor field A on TM \ {0}, its dynamical covariant derivative is given by 

(2.6) VA = VoA-AoV. 

From first formulae in (|2.4p and (|2.6p it follows that Vh — and Vv = 0. Hence V preserves the 
horizontal and vertical distributions H and V. However, V does not act in a similar way on these 
two distributions. This can be seen locally as follows. Using formulae ()2.4p and ()2.5p we have 

(2.7) = A^f Vdx^ = -N^dx', 



dt \ dt 'Jdyi' ^ \ dyi 
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The action of the dynamical covariant derivative V on the components of a horizontal vector field 
X ~ hX is given by 

ox'- J ox^ 



(2.8) VX = V(X'— ) =VX^ — , VX^^SiXn + N^X^. 



In formula (|2.8p both components X' and VX* transform, under a change of induced coordinates 
on TM, as the components of a vector field from the base manifold A/. The first order differential 
operator VX' was introduced by Kosambi [27] with the name of bi-derivative. The term "dynamical 
covariant derivative" was introduced by Carihena and Martinez in |14) . 

Next lemma gives some compatibility conditions between some of the geometric structures 
introduced so far. Let us emphasize that in the next lemma we do not yet assume any relation 
between the semispray S and the nonlinear connection {h, v). However, we will use this lemma to 
fix a nonlinear connection in Proposition 12.71 

Lemma 2.1. Consider a semispray S, a nonlinear connection {h,v), and the dynamical covariant 
derivative V associated to the pair (S, {h,v)). Then, 

(2.9) h o Cg o J = —h, J o Cs o V = —v, 

(2.10) CsJ + ld-2v = VJ, \7J = 2(—--Nn—(g,dx\ 

V dy' J dy3 

Proof. For an arbitrary vector field X G X{TM \ {0}) we have X + [S,JX] & X"(TM \ {0}). 
Composing from the left with h it follows that hX = h[JX,S], which shows that first formula 
(|2.9p is true. Similarly, composing by J, setting JX = vZ for some Z G X{TM \ {0}) gives the 
second equality in formula (|2.9p . 

Using formula (|2.4[) and formulae (|2.9p we obtain that 

V o J = w o Cs ovoj^vo Cs o J = (Id —h) o Cs o J ~ Cs o J + h, 
JoV = JohoCs°h = Jo Cs o h ^ J o Cs o (Id —v) = J o Cs + v. 



Using the above two formulae and formula (|2.6p we obtain that VJ = CsJ + Id— 2u, which 
is formula (|2.10p . Second formula (|2.10p follows directly from the local expressions (|2.7p of the 
dynamical covariant derivative. □ 

We note that from formula p.lOp . the (l,l)-type tensor VJ = CsJ + Id— 2t; represents the 
difference between an arbitrary nonlinear connection {h,v), which is not fixed yet, and the canonical 
nonlinear connection associated to a semispray S that will be fixed in Proposition 12.71 

2.2. Symmetries for systems of SODE. It is well known that a semispray induces a canonical 
nonlinear connection which in turn, determines the five geometric invariants of the semispray. This 
is known as KCC theory, after Kosambi, Cartan, and Chern |2]. In this section we show that this 
connection can be fixed using symmetries of the semispray. This approach follows the statement of 
Kosambi [28] "the connection cannot be determined without recourse to the equations of variation". 
In the next section we show that this approach of fixing a connection also generalize to systems of 
higher order differential equations. 

Definition 2.2. A vector field X 6 X{TM \ {0}) is a dynamical symmetry of a semispray S if 
[S,X] =0. 

A direct calculation shows that a vector field on TM \ {0}, locally expressed as 
(2.11) x = X\x,y)^^+Y^ix,y)^^, 
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is a dynamical symmetry if and only if Y"^ = and 
(2.12) S'\X')+X{2G')^0. 

It follows that, while studying dynamical symmetries for a semispray S, the following set of vector 
fields on TM \ {0} plays an important role: 



(2.13) 



Xs-[xe X{TM \ {0}), X = A + 5(X^) A| 



A vector field X e is called a newtonoid [35]. Without local coordinates, the set of newtonoid 
vector fields can be expressed as follows 

(2.14) Xl = Kcr (J o Cs) = Im (Id+J o Cs) ■ 

For a vector field X = X^d/dx^ G X{M) the complete lift, X^-^ G X{TAI \ {0}), is the vector 
field defined by 

r) BX^ r) 

(2.15) ^"=^'(^)^+a^(^)^V- 

The vertical lift of a vector field X G X{M) is the vector field X^'^ G X{TM \ {0}) given by 
Xi.o = J(Xi'i) gg. We denote by X^'^{TM \ {0}), the set of complete lifts of vector fields on M, 
and by X^^°{TM \ {0}) = J {X^'^{TM \ {0})), the set of vertical hfts of vector fields on M. 

From expression p.lSp and formula (|2.15p it follows that complete lifts and newtonoid vector 
fields are related by 

X^^^{TM\{0})C Pi X\, 

S scmispray 

and equality holds if and only if dim A/ > 2. If dim A/ — 1 the above inclusion is strict due to the 
fact that the fibers of TM \ {0} are not connected. 

Definition 2.3. A vector field X G X{M) is a Lie symmetry of a semispray S if its complete lift, 
X^'^ , is a dynamical symmetry, which means that [S", X^'^] = 0. 

Next, we define a C°°-module structure on the set of newtonoids. This module structure was 
introduced in |40) for discussing adjoint symmetries for systems of ordinary differential equations. 
In proposition 12 . 71 we show that this structure can be used to characterize the canonical nonlinear 
connection associated to a semispray. 

Remark 2.4. For / G C°°{TM \ {0}) and X G X^, we define the product 

/ * X = (Id + J o Cs) ifX) = fX + S{f)JX. 

i) With respect to the product *, the set Xg has the structure of a C°°{TAI \ {0})-module, 
and the set X'^^'^{TM \ {0}) has the structure of a C°°(M)-module. 

ii) The maps J : (X^,*) (X''(TM \ {0}),-) and h : (X^,*) -> (X''(TA/ \ {0}),-) are 
isomorphisms between C°°{TM\ {0})-modules. The map Dn : (X^'i(TAf \ {0}), *) ->■ 
(X(Af), •) is an isomorphism between C°°(A'/)-modules. 

iii) A vector field X on TM \ {0} is a newtonoid if and only if it can be expressed as follows 

X^X^i.,y).^. 

We have seen that a vector field X on TM \ {0} is a dynamical symmetry if and only if it is a 
newtonoid and satisfies equation (|2.12p . Therefore, a vector field X = X'^{x)d/dx^ on M is a Lie 
symmetry if and only if components X^ satisfy the system of equations (j2.12p . Our aim now is 
to rewrite equations (|2.12p such that its terms will have a covariant character. Note that neither 
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S'^(X'), nor X(2G'*) in formula (|2.12p . has such a covariant character. For this we evaluate first the 
horizontal and vertical components of a vector field, which behave as the components of a vector 
field from the base manifold. This will allow us to characterize newtonoid vector fields using the 
dynamical covariant derivative. 

Lemma 2.5. Consider a semispray S, a nonlinear connection {h,v), and the dynamical covariant 
derivative V associated to the pair (S*, {h, v)). A vector field X on TM \ {0} is a newtonoid if and 
only if v{X) = J(VX), which locally is equivalent to 

(2.16) ^ = + 

for some functions defined on the domain of induced coordinates on TM \ {0}, where VX* is 
defined in formula (j2. 



Proof Since J oV = J o Cs + v it follows that J[S,X] = if and only if v{X) = J{VX). The 
local formula ((TTC)) follows by formulae (gH) and ^J^. □ 

Proposition 2.6. Consider a semispray S, a nonlinear connection {h,v), and the dynamical 
covariant derivative V associated to the pair {S, {h^v)). A vector field X £ X{TM \ {0}) is a 
dynamical symmetry if and only if X is a newtonoid and satisfies 

(2.17) V(JVX) + <I>(X) = 0. 

Proof. A vector field X e X{TM \ {0}) is a dynamical symmetry if and only if h[S,X] ~ and 
v[S,X] = 0. Composing by J shows that the first condition is equivalent to J[S,X] = 0, which in 
turns is equivalent to X £ Xg. If X £ X{TAI \ {0}) we have 

v[S, X] = v[S, vX] + v[S, hX] = y{vX) + <P{X). 

Thus X is a dynamical symmetry if and only if X G Xg and \7{vX) + ^{X) ~ and the result 



follows by Lemma [2751 □ 

We note that in the first term of formula (|2.17p . JVX only contain horizontal derivatives of 
X, while the outer V acts on the vertical vector field JVX. According to formula (|2.7p . the two 
actions might be different. The next proposition shows that the two actions of the dynamical 
covariant derivative V coincide if and only if V J = is satisfied. Moreover, this condition is only 
satisfied for one nonlinear connection. In terms of S this unique nonlinear connection is explicitly 
given through its vertical projector 2v = Id+£sJ, which is the canonical nonlinear connection in 
|24| . The advantage of this method of assigning a nonlinear connection to a semispray is that the 
method also generalize to systems of higher order ODEs. This is the topic of Section |31 

Proposition 2.7. Consider a semispray S, a nonlinear connection {h,v), and the dynamical co- 
variant derivative V associated to the pair (5*, {h, v)). Then, the following conditions are equivalent: 

i) V restricts to a map V : Xg Xg that satisfies the Leibniz rule with respect to the * 
product; 

ii) V J = 0; 

iii) CsJ + ld-'2v ==0; 

iv) AH = dG'/dyK 

Proof. The formulae in (j2.10p show that conditions ii), iii) and iv) are equivalent. 

To prove implication ii) =► i), let X £ X\. According to Lemma [^75l we have that vX = JVX. 
We apply V to both sides of this equality, use Vw = and V J = 0, and obtain v (VX) = VuX = 
V (JVX) = JV (VX) , which according to Lemma shows that VX G X\. 
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Let / G C°°{TM \ {0}) and X e X^. Since f*X^fX + Sif)JX and V satisfies the Leibniz 
rule for the • product, it follows that 

(2.18) V(/ * X) = S{f)X + fVX + S\f)JX + S{.f)V{JX), 
V(/) * X + / * VX = S{f)X + fWX + S^{f)JX + S{f)J{VX). 

Using the condition V J = we obtain V(/ * X) ^ V(/) * X + / * VX . 

To prove implication i) =^ ii), we will show that the (1, l)-type tensor field VJ vanishes on 
the set X\ U X^iTM \ {0}), which is a set of generators for X{TM \ {0}). By formulae ((^ 
and (|2.6p . VJ vanishes on X^(rA'/ \ {0}). Since V satisfies the Leibniz rule with respect to the 
* product, it follows from the two formulae (|2.18p that S{f){V J)X = for an arbitrary function 
/ e C°°{TM \ {0}) and arbitrary vector field X G X^. Therefore, VJ = on X^, and this 
completes the proof of the proposition. □ 

For a semispray S, we say that the nonlinear connection characterized by any of the four 
conditions of Proposition 12.71 is the canonical nonlinear connection of S. For the remaining of 
this section we will consider only this canonical nonlinear connection induced by a semispray S. 
The next proposition characterizes the dynamical and Lie symmetries of a semispray S using the 
dynamical covariant derivative V associated to S and its canonical nonlinear connection. 

Proposition 2.8. Consider a semispray S . 

i) A vector field X S X(TM \ {0}) is a dynamical symmetry of S if and only if it is a 
newtonoid and 

(2.19) VVX + 4'(X) = 0. 

ii) Locally, a newtonoid vector field X G X(TAI \ {0}) is a dynamical symmetry of S if and 
only if its horizontal components satisfy 

(2.20) V'^X' + R]X^ = 0, 

wtiere VX* is defined in equation (|2.8p . 

iii) A vector field X G X(A/) is a Lie symmetry of S if and only if 

(2.21) V'^X^° + ^{X^^^) = 0. 

Proof. Since we are using the canonical nonlinear connection, we have that VJ = 0. The first 
claim follows using Proposition 12.61 For the second claim we use local formulae (|2.8p and (|2.3[) . 
Therefore, the local components of the vertical vector field W"^ JX + ^{X) are V^X* + i?* X-' . Hence 
formulae p.l9p and (|2.20|1 are equivalent. For the third claim consider a vector field X G X(M). 
Its complete lift X^-^ is a newtonoid for S and formula (|2.2ip follows from formula (|2.19p since 
JXi-i=Xi'0. □ 

Equation (|2.20p represents the covariant reformulation of equation (|2.12p . both terms V^X* and 
KjX^ behave as components of a vector field from the base manifold. 

2.3. First order variation and Jacobi fields for systems of SODE. Consider c : / — > AJ 
a geodesic of a semispray S. A geodesic variation of c is a smooth map V : I x (— £, e) — > M, 
V = V{t,s) such that 

i) V{t,0) c{t), for all t in /, 

ii) V{t, s) is a geodesic of S for all s in (— e, e). 

Definition 2.9. A vector field ^ : / ^ TM, along a geodesic c = tt o ^ of a semispray S, is called 
a Jacobi field of S* if ^ is the variation vector field, ^(t) — dsV(t, 5)1^^0, of some geodesic variation 
V. 
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Next we characterize the Jacobi fields of a semispray, and hence the geodesic variation, in terms 
of a system of differential equations, using the dynamical covariant derivative and the curvature 
components of the Jacobi endomorphism. We define first the dynamical covariant derivative of a 
vector field along a curve. 

Consider (, : I TM, ^{t) = {x' (t) , C (t)) , a vector field along a regular curve c ~ ir o ^. Then 
the complete lift of £, is the vector field ^^-^ : / -> T(TM \ {0}) along c', defined as ^^'^(0 = 
{x'{t), dx'/dt, ^'(t), dC/dt) . The vertical lift of ^ is the vector field C^-O : / T{TM \ {0}) along 
c', defined as = [x''{t),dx^ /dt,Q,^^{t)y Since ^ is a vector field along c, for each t ^ I, 

we can find a vector field X S X{M) such that X o c = ^ near t. Then ^^'^ = X^'^ o c' and 
= o c'. Assume now that c is a geodesic, which means that c' is a geodesic of a semispray 
S. Then, we define 



o c 



(2-22) Vf = + A^; ( X, — ) e 



If locally ^ = {x\C), then V^^'" = VCd/dy\ where 

dt ^ \ dt. 
Thus V^^^° does not depend on the choice of X. 

Proposition 2.10. Consider a bounded interval I. A vector field ^ : / — TM, along a geodesic 
c = TT o ^, is a Jacobi field if and only if 

(2.23) V2^i'0 + $(^i'i) =0. 

If locally ^{t) = (x*(i), ^'(t)), equation (|2.23p is equivalent to 

(2.24) vY + i?}(^W,§)e^=0. 

Proof. Let ^ : / — > TAf, ^(i) = (a;*(t), ^'(t)) be a vector field along the geodesic c = tt o ^. Using 
formula (|2.22p . we have 

(2.25) V2^i'° + $(e''i) = 

If ^ is a Jacobi field then there is a geodesic variation V{t, s) such that V{t, 0) = c{t), 
— (i,0) = f (i), and ^ + 2G' f — j = 0,V(i,s) G / x (-£,£). 

Differentiating the last equation with respect to s and setting s — we obtain that the right hand 
side in formula p.25p vanish and equations p.23p and (|2.24p follow. 

Conversely, if ^ : / — TM, ^(t) = (a;*(t), ^'(<)) is a solution of the Jacobi equations (|2.24|) then 
right hand side in formula (|2.25p vanish. Using the assumption that the interval / is bounded and 
[TT| Theorem 4.4] it follows that ^ is a variation vector field for some geodesic variation. Therefore 
^ is a Jacobi field. □ 

Consider c{t) = {x'^{t)) a geodesic of a semispray S. Its tangent vector field c'(t) = {x'^{t), dx^ /dt) 
is a Jacobi field of S. This can be seen either from the fact that the right hand side in formula 
(|2.25[) vanish for = dx^ /dt or from the fact that c' is the variation vector field of the geodesic 
variation U(t, s) = c{t + s). In the Riemannian or the Finslerian context [4j [15] we can also see 
that the tangent vector field c' of a geodesic c satisfies the Jacobi equations p.24p and hence it is a 
Jacobi field since the geodesic equations <\2.2\i can be written as V {dx^ /dt) ~ and the curvature 
components satisfy i?^- (a;, dx/dt) dx^ /dt = 0, [421143], 



For a geodesic c{t) ~ (a;*(t)) of a semispray S, the vector field {x^{t),tdx^ /dt) is a geodesic of S 
if and only if y^dC^/dy^ ~ 2G". This condition means that functions G* are 2-homogeneous in the 
fibre coordinates y. In this case S is called a spray and this includes Riemannian and Finslerian 
cases. 

The second implication of Proposition !?. lOl is proved in the Riemannian |15) or Finslerian context 
[4] using the exponential map. In the general context of a semispray. due to the lack of homogeneity, 
this exponential map cannot be defined and the proof in jlll Theorem 4.4] uses the fiow of the 
complete lift of a semispray. 

In the covariant form of the Jacobi equations (|2.24p . the components i?* of the Jacobi endo- 
morphism contain information about the geodesic behavior and the stability of the geodesies of a 
semispray S. If the real part of the eigenvalues of the curvature components i?* are negative then 
the geodesies will spread out and will be unstable, while if they are positive the geodesies will be 
stable ISg. 

3. Symmetries and first order variation for systems of higher differential 

equations 

In this section we will extend the geometric theory for systems of SODE, developed in the 
previous sections, to systems of HODE. A system of (fc + 1) order ordinary differential equations 
can be associated with a semispray S of order k, which is a special vector field on the fc-th 
order tangent bundle. For a semispray of order k there are various nonlinear connections one can 
associate to it [U El [121 IHl IHl EQ] . In the first part of the section we consider a semispray of 
order k and an arbitrary nonlinear connection. For this pair we define the dynamical covariant 
derivatives of first and higher order, and the Jacobi endomorphism. The main contribution is to 
describe a particular nonlinear connection that yields covariant expressions for symmetries and 
first order variation of the system of HODE. For this nonlinear connection we explicitly compute 
all components of the Jacobi endomorphism and use these components to provide a geometric 
expression for the Wuenschmann invariant [20l [2TJ [38] . 

3.1. Systems of HODE and geometric structures on T^Ad. The framework for studying 
systems of (fc + 1) order ordinary differential equations on a manifold M is the tangent bundle 
T'^M = J^M of order > 1 [HI |30l |34l Sll • This is the fc-th order jet bundle of curves c from a 
neighborhood of e K to M. For a curve c: I ^ M, c{t) = {x\t)), denote by j'^c : / T^M, its 
fc-th jet lift, defined as 

I dx\ , 1 d^x^ 



If c is a regular curve, then j^c{t) £ T^M \ {0} for all t G I. Local coordinates on M are denoted 
by (a;') and the induced local coordinates on T'^M are denoted by {x\y'^'^^\ • ■ • , y^^')'), where 



y ^^""^^-a! dt<- 



a G {l,..,fc}. 



Let also Denote T^M ~ M . For each a e {0, 1, k — 1}, the canonical submersion 

TT% : T^M T"M induces a natural foliation of T'^M. Consider also the subbundle T'^M \ {0} = 
&T''M,yW ^0}. Note that T'=M \ {0} = (ttJ^)"^ (TAf \ {0}). 
The tangent structure (or vertical endomorphism) of order k is the (1, l)-type tensor field on 
T'^M defined as 

It has the following properties: = 0, Im J" = Kcr a e {1, k}. 
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The foliated structure of T^M determines k regular vertical distributions 

Vdu) = KcrDuT^l^i = Im = Ker J^-"+\ for u G THl^aG {1, k}. 

For each a G {1, fc}, the distribution is tangent to the fibers of tt^_^i : (x*, y'-"'^'*, • ■ • , yC^')*) 
(cc', y^^-**, • • • , y'-""^-''), and hence it is integrable. Note that diml/o. = (fc — a + l)n, a G {1, k} 
and Vfe(u) C Vfc_i(M) C • • • C for each u e T'^M. 

In this section we pay attention to the geometry of a system of (fc + 1) order ordinary differential 
equations. Such a system of HODE on M can be represented using a semispray of order k, which 
is a special vector field on T'^M \ {0}. 

A semispray of order fc on Af is a vector field S G X{T^M \ {0}) such that any integral curve 
of S", 7 : / T^M \ {0} is of the form 7 ~ 3^{t^q ° 7). In induced coordinates for T^M \ {0}, a 
semispray of order k is given by 

(3.2) 5 = y^^^^7^ + 2y(^^'7rL + --- + fcy^'=^'^-7^-(fc + l)G^ ' 



for some functions G" defined on domains of induced local charts. For an integral curve 7 ; / ^ 
T^AI \ {0} of S', we say that curve c = ttq o 7 is a geodesic of S . Therefore, a regular curve 
c : / — > M is a geodesic of S if and only if S* o j^c — {j^c)' . Locally, a regular curve c : I M, 
c{t) = {x^{t)), is a geodesic of 5* if and only if it satisfies the system of {k + 1) order ordinary 
differential equations 

1 d^^^x^ f dx 1 d'^x\ 

^^■^^ (fc + 1)! dffc+i + ' ■ ' ■ ' fc! dt^ j " 

Thus a semispray of order k describes systems of HODE with regular curves on M as solutions. 

Following the line developed in Section [2] we will show that the geometry of a semispray of 
order k can be determined from a canonical nonlinear connection, associated to the semispray. An 
n-dimensional distribution ifo : u G T^M \ {0} ^ Hq{u) C TuT^M, that is supplementary to the 
vertical distribution Vi{u), is called a nonlinear connection, or horizontal distribution. Therefore, 
for a nonlinear connection, we have the following decomposition TuT^M = Hq{u) © Vi{u), for 
all u G T'^AI \ {0}. For a nonlinear connection Hq, consider also the n-dimensional distributions 
Ha = J"{Hq), forQ:G{l,...,fc — 1}. Therefore, for a nonlinear connection, we have the following 
decomposition: 

(3.4) TuT'^M = Q'^.ZlHaiu) ® Vkiu),yu G T'^M \ {0}. 

The set of k distributions Hq, Hi, Hk-i is called a multi-connection in |4l]. We will denote by 
ho,hi, ...,hk-i,Vk the projectors that correspond to decomposition p.4p . which characterizes the 
nonlinear connection. Sometimes, for uniform notations, we also denote hk = Vk- Decomposition 
induces also a decomposition of the vertical subspaces 



(3.5) Vaiu) = ©yjaHpin) © Vk{u),'iu G T'^M \ {0}, Va G {1, k - 1}. 

Decomposition (|3.5p implies that Idjimjo — Id iKorjfc-o+i ~ X]^=q^/3i ^ {0,...,fc}. 

Therefore, for all a G {0, fc} we have the following identities 

k a — 1 

(3.6) J" = ^/i^oJ", J^-"+i = ^ J^-^+Iq/i^. 

I3=a /3=0 

For a, /3 G {0, fc}, one can reformulate the above identities as follows 

0, if^<a-l, 
J" o hfj-a, if /? > a, 
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(3.7) /i^ o J" = 



(3.8) rohp^ 



0, ifa + /3>fc + l, 

/ia+/3oJ", ifa + /3<fc. 

For a e {0,...,fc - 1}, we denote by X'*° (T'^A/ \ {0}) = K (X(r'=M \ {0})), X^*- (T'^A'/ \ {0}) = 
Vk {X{T^M \ {0})), the C°°{T''M \ {0})-niodules of horizontal and respectively vertical vector 
fields. 

For a semispray of order fc, S, and a nonlinear connection (ha, Vk), consider the vertically valued, 
(1, l)-type tensor field on T'^M \ {0} 

fc-1 fc-1 
(3.9) $ = -Vk O CsVk = Wfc o ^ £s/ia = Vk o ^ (£s ° ha) 

a=0 a=0 

that we will call the Jacobi endomorphism. Here, for a (1, l)-type tensor field A on T^M \ {0}, 
= Cso A — Ao Cs denotes the Frolicher-Nijenhuis bracket of S and A. We will obtain now the 
local expression of the Jacobi endomorphism. On the domain of induced local charts on T'^A/\{0}, 
consider the vector fields 

for a G {!,..., fc}, which form a basis for ?i{T^M \ {0}), adapted to the decomposition p.4p . 
Functions -/V^'q,)j are locally defined on T'^Af \ {0} and are called the coefficients of the nonlinear 
connection. The dual basis to the basis p.lOp is given by the following locally defined 1-forms on 
T^M\{G}: 

a 

(3.11) dx\ 6y^'^^' = dy^'^^^ + ^ Af('^)/y("-«^ a G {1, k}, 

13=1 

where the dual coefficients M'l^p-^j are given by |34| 

Q-l 

(3.12) Af(\)^. = Nl,^^, Ml^^ = Nl^^ + Y: NI-P)sMIp),, a G {2, k}. 

A change of induced local coordinates (a;',y'^"'*) — >• (p(x),y^"^* (x,?/*^^-')) on T^M induces the 
following transformation rule for the basis p.lOp and its dual basis p. lip 

5 9P 6 5 dx^ 5 d dx^ d fl fc 1) 

= —dx\ Sy^'^^' = d^^y^''^'' a G {1, 

Therefore the components of a tensor field on T'^AI, with respect to these bases, will transform as 
the components of a tensor field on the base manifold M. 

The projectors {ha,Vk) that correspond to the decomposition p.4p can be expressed, with 
respect to the bases p.lOp and p. lip as follows 

5 S d 

= ^ ® dx\ ha = ^ ® a G {1, k - 1}, Vk = ^ ® <5yW^ 

With respect to the bases p.lOp and p. lip the tangent structure of order k has the following 
expression: 

(3.13) J = — ^ ® dx' + 5y^^^' + ■■■ + —Ty ® 5y^''-^^\ 
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Locally, the Jacobi endomorphism can be expressed as follows 



(3.14) 



a=0 



® Sy 



Under a change of induced local coordinates on T^M^ the k components R\^g,^j of the Jacobi 
endomorphism transform as the components of a (1, l)-type tensor field on the base manifold. 
We will refer to R\a)j '"^^ curvature components. In Subsection 14.11 we will show that in 

the Riemannian context, these components are functions of the Riemannian curvature and its 
dynamical covariant derivatives. 

We will provide now explicit formulae for the curvature components R}^a)j of the Jacobi endo- 
morphism. For a semispray of order k, S, and a nonlinear connection with coefficients N^^^j we 
introduce the following notations: 



a-l 



(3-15) = aA^(^„),-5(iV('_i),.) + ^(i),^(.-i). + E^(«.^(-«- 

/3=2 

for a G {3,...,fc}. With respect to these notations we have the following formulae for the Lie 
bracket of the semispray S and the vectors of the basis p.lOp 



5x^ 



(3.16) 



'6x' ' ^ msyi0-^)^ 

13=2 " 

5G 



d 8 k — a 

------- + ^(ib-:w^ + E^( 

B=2 



13=2 
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{13)3 gy{a+P~l)i 



SG 



E ^03)i^(fc+l-a- 



Qy(k)i ' 



for each a e {1, k - 1}. From formulae ((XTC)) it follows that 6y'^°'^'[S, 5/6y'^^^^] behave as (1, 1)- 
type tensors on M, for all a ^ G {0,1, ...,fc}. Therefore, ^ {2,---,fc}, behave as 

(1, l)-type tensors on A/. 

From formulae (|3.9p and (|3.14p it follows that the curvature components of the Jacobi endo- 
morphism can be computed as follows 



R 



d 

Qy{k)^ 



Vk 



e {o,i,...,fc-i}. 
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Using the expressions p. 161) we obtain the foUowing formulae for R^i^^yj-, for a 6 {0, 1 



1} 



R) 



SG 



k — a 



(3.17) 



(/3)j^(\+l-a-^)s- 



3=2 



For a fc-semispray S and a nonUnear connection (/iq,, w^), consider the map V : X(r''Af \ {0}) 
X(r'=A/\{0}), given by 

fe-i fc-i 
(3.18) V = ^ /la o £5 o /ic, + ufc o £5 o £5 + ^ /iq, o + Vk o CsVk- 

a=0 a=0 

We call V the dynamical covariant derivative associated to the pair (S*, (/iq, u^-)). By setting 
V/ = S{f), for / G C°°{T''M\{0}), using the Leibniz rule, and the requirement that V commutes 
with tensor contraction, we extend the action of V to arbitrary tensor fields and forms on r''7\/\{0}. 
Formula (|3.18p implies that V/iq, = for every a € {0, fc — 1}. Hence V preserves all distributions 
Ha. Similarly, we obtain that Vufc = 0, hence it preserves the vertical distribution Vk- From 
formula (|3.5p it follows that V preserves also all the vertical distributions Va ■ Dynamical covariant 
derivative V acts in a similar way on the distributions Ha, but has a different action on the vertical 
distribution Vk- This can be seen locally as follows. Using formulae p.l6p and (|3.18p we obtain 



(3.19) 



d 



Vk 



d • 



a e {0,...,fc- 1} 



= (fc + 1) 



(1)1 



d 



The action of the dynamical covariant derivative V on a horizontal vector field X = haX, for 
a S {0, 1, k — 1}, is given by 



(3.20) 



by 



yx' 



x^. 



The first order differential operator VX* and its higher order iterations were considered by Kosambi 
|28j . For k = 2, iterated actions of V where also considered in [5] in order to fix a nonlinear 
connection for a second order Lagrange space. 

Next lemma gives some compatibility conditions between the geometric structures introduced 
so far. In this lemma we do not assume any relation between the semispray and the nonlinear 
connection. However, we will use this lemma to fix partially the nonlinear connection later on. 

Lemma 3.1. Consider a semispray S, a nonlinear connection {ha,Vk), and the dynamical covari- 
ant derivative V associated to the pair (S, (ha, Vk))- Consider a, /3 G {0, k}. Then, the following 
formulae are true: 



(3.21) 
(3.22) 

(3.23) 
(3.24) 



O Cs O ha 



'{13 + l)hpo Jf^, ifP^a-l 
0, ifP<a-l, 



-aJ^-" o ha, 
0, 



ifa = k-l3 + l 
if a > k - 13 + 1, 



V J" o J'' 



k-c 



v.r = LsJ 



k 



{khk-1 
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Wfc)o J^-\ i/a> 1, 



(3.25) vj" = (fc + l)(^^-iV(',)^.j^0<Sy('=-"^ ^fa>l, 

(3.26) CsJ + ld-{k + l)vk = 

Proof. From formula (|3.ip . a local computation shows that for every a G {1, k} and for every 
X e X{T''M \ {0}) we have 

(3.27) [S, J"X] - J^iS*, X] + ar-^X e Kcr J = Im 

For a G {!,..., fc}, formula p.2ip follows by composing in formula (|3.27p to the left with J. For 
a = 0, formula (|3.2ip follows since £s Id = 0. 

Formulae (|3.22p follow from formula (|3.27p by composing to the left with hp and using the 
identities ([3T7|) . 

Formulae (|3.23p follow from formula (|3.27p by composing to the right with hj3 and using the 
identities p.Sp . 

From formulae (|XT5| . (jXT)) . (j^ . (jX^ . and ((X^ for ;3 = /c, we obtain 

VoJ*^ = VkoCsoj'^^CsoJ^-{ha + --- + hk-i)oCsoj'^ 

= Cs ° J"" + khk-1 o 

J*-^ O V J*-' O Cs O ho ^ j'' O Cg - j'' O Cg O {hi -\ h + Wfc) 

Now the second equality in formula p.24p follows using the fact that VJ'^ = V o - oV. 

We wiU prove now the first equality of formula (|3.24p . Since VJ"o j'^^" = Vo j'^ — J"oVo J*^"", 
it remains to show that J" o V o J*' o V, for a G {1, fc}. Using formula p.lSp we have 

J" o V o j'^-" = J" o hk-a oCgO hk-a o j''"" = Wfc o J" o Cg o o /iq. 

In formula p.2ip . we replace a by fc — a, compose to the left with J"^^ and obtain 

roCgo J'^-" = J^oCg- {k - a)J^-^. 

Using the above two formulae and the fact that Vk o J^^^ o /ig = 0, we obtain 

J" o V o J*=-" = oCgo ho = o V. 

Formula (|3.25p follows from expression (|3.13p and the local expression of the dynamical covariant 
derivative (|3.19p . 

By direct calculation it follows that for a G {0, fc — 1} 



Using these formulae and p.25p for a 1 it follows that formula (|3.26p is true. □ 

When = 1 formulae ((X^ and read as follows 

(3.28) ho o Cg o J = —ho, J o Cg o Vk ~ —kvk- 

Formula (|3.26p generalizes formula (|2.10p . In the case k = 1 the tensor VJ uniquely determines the 
tensor Z^sJ + Id— 2w that fixes the canonical nonlinear connection. For the case fc > 1, a nonlinear 
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connection can be fixed using the (l,l)-type tensor CsJ + Id —{k + l)ufe. Using formula p.26p . in 
Theorem 13. 101 we will determine this tensor by fixing VJ and 



Therefore, in Theorem 13. 10[ we will determine a nonlinear connection for a semispray S by fixing 
VJ and the (1, l)-type tensor whose local expression is given in formula p.29p . 

3.2. Symmetries for systems of HODE. For a semispray of order k there are various nonlinear 
connections one can associate to it [1] |6l [161 IHl ISS] . In this section we will show what are the 
advantages of using the nonlinear connection proposed by Miron and Atanasiu [36] . We show that 
this connection is very useful for studying the symmetries and the first order variation of a system of 
higher order ordinary differential equations. Moreover, we show that with respect to this nonlinear 
connection one can explicitly write down all components of the Jacobi endomorphism. We use 
these components to provide a geometric expression for the Wuenschmann invariant |20[ [2T1 [38] as 
well as for an invariant introduced by Fels |22| for the inverse problem of a fourth order ODE. 

Definition 3.2. A vector field X e X[T^M \ {0}) is a dynamical symmetry of a /c-semispray S if 
[S, X] = 0. 

A direct calculation shows that a vector field on T'^M \ {0}, 

(3-30) x^x^ ^ + 1 ,(^)) 

a— 1 ^ 

is a dynamical symmetry if and only if its components satisfy alF*-"-** = S"{X''), for all a € 
{1, k} and 

(3.31) S''+^X') + {k + l)\X{G')^0. 

It follows that, while studying dynamical symmetries for a fc-semispray S, the following set of 
vector fields on T'^'A/ \ {0} plays an important role: 

(3.32) X| = |x e XiTHj \ {0}),X = X' A + . 

A vector field X S is called a newtonoid of order k. Next, we will provide global characteri- 
zations, in terms of some differential operators for the set of newtonoid vector fields X| and their 
images J" for a e {1, fc}. 

Consider 7r| : X{T^M \ {0}) — > X{T^M \ {0}) the differential operator of order k defined as 



's 

4 = Id +i / o £5 + . . . + 1 J^- o 4. 



The form of the differential operator tt^ was inspired by the generalized Cartan operator acting on 
1-forms on T^M, considered in [TBI I40| . 

Lemma 3.3. For each a G {0, the set J" can he expressed, without local coordinates, 

as follows 

(3.33) r {X%) = Ker ( J o £5 + aid) n Im J" = Im(J" o ti%). 
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Proof. First, we prove the two equalities in formula (|3.33p for a = 0. Consider a vector field X on 
r'"'A/\ {0}, given locally by formula (|3.30p . It follows that J[S, X] ~ ii and only if its components 
satisfy /3!r(^)' = S^{X'), for aU /3 e {1, k}. This proves first equality in formula (p33| . 

From formula (j3.2ip we have J o Cs° J = o Cs — J ■ Using this formula recurrently, we obtain 
that (J o Cs) o 7r| = 0. Therefore, Im7r| C Ker (J o Cg). 

For the converse inclusion, consider a vector field X on T'^M \ {0} such that J[S', X] ~ 0. It 

follows that C^X G Ker J/^, for all /3 e {1, k}. Hence, ( J'^ o X = and X = 7r|(X). 

Consider now a G {1, fc}. Using the above considerations it follows that J" (-^1) = Im(J" o 
7r|). We prove now the first equality in formula (|3.33p . Consider Z = J°'X for X G X|. Therefore 
J[S,X] = and using formula we have J[S', J"X] + aJ"X = J°'+^[S,X] = 0. It follows 

that Z G Ker(Jo£5 + ald) nim J". 

For the converse inclusion, consider Z = J^X for X G X(T'=M\{0}) such that J[5, Z] +aZ = 0. 
Using again formula (|3.2ip it follows that J"+-'^[S', X] =0. li X is given locally by formula 
(fOO]) . then its components satisfy /Jiyt-^)* = ^'^(X'), for all /3 G {l,...,fc - a}. It follows that 
Z = J"7r|XG J"(X|). □ 

For a vector field X ^ X^d/dx' G X(M), the complete lift, X^^^ G X{T^M \ {0}), is the vector 
field defined as 

(3.34) x''^'' ^X\x)^ + y -d'^iX'ix))-^, 



where dx is the Tulczyjev operator 



For each a G {0, fc — 1}, the vertical lift of a vector field X G X(M) is the vector field X*^'" = 
jk-a (^x^,k^ G X(r'^7\/ \ {0}), [Ullin]- We denote by X^^^{T^M \ {0}), the set of complete lifts 
of vector fields on M and by X'='"(r'=M \ {0}), for a G {0, k - 1}, the set of vertical lifts of 
vector fields on M. From expressions (|3.32p and p.34p it follows that complete lifts, vertical lifts, 
and newtonoid vector fields are related by 

(3.35) X'=^"(r'^'Af \{0}) C Pi J'^--"(X^), 

S semispray 

for all a G {0, fc}, and the equality holds if and only if dim Af > 2 and a = k. If dimi\/ = 1 the 
above inclusion is strict due to the fact that the fibers of T'^M \ {0} are not connected. 

Definition 3.4. A vector field X G X(M) is a Lie symmetry., or a point symmetry, of a fc-semispray 
S if the complete lift X^'^ is a dynamical symmetry, which means [S*, X'^'-'^'J ~ 0. 

Next, we define a C°°-module structure on the set of newtonoid vector fields of order k. This 
structure is inspired by a module structure on a set of 1-forms associated to adjoint symmetries of 
a system of HODE, [ID]. In proposition 13.81 we show that this structure can be used to fix a part 
of a canonical nonlinear connection one can associate to a semispray. 

Remark 3.5. For / G C^{T^M \ {0}) and X G X|, we define 

k 

(3.36) /*X = 7r|(/X) = V-5"(/)J"X 

— ' a\ 

a=Q 
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i) With respect to the product the set J" has the structure of a C°°{T^M \ {0})- 

module and the set X^'^^iT^M \ {0}) has the structure of a C°°(A/)-module, for all a £ 
{0,1, ...,fc}. 

h) The map K ■ ( J" , *) ^ (X''° {T'^M \ {0}), •) is an isomorphism of C°°{T''M \ {0})- 
modules, for all a S {0,...,fc}. The map Dtt^ : {X''-''{T''M \ {0}),*) ^ (X(M),-) is an 
isomorphism of C°°(M)-modules. 

iii) A vector field X on T'^M\{0} belongs to the set J" (X^) if and only if it can be expressed 



{^'y'''--y^'')*d^^ vae{o,i,...,fc}. 



As we have seen, a vector field X on T'^AI \ {0} is a dynamical symmetry if and only if it is 
a newtonoid and satisfies equations p.3ip . Our aim now is to rewrite equations (|3.3ip such that 
its terms have a covariant character. Note that neither one of the two terms of equation p.3ip . 
gk+if^j^z-^ and X(G"), has such a covariant character. To accomplish this goal, we first characterize 
newtonoid vector fields in terms of their expressions in the adapted basis p.lOp . The components 
of a newtonoid vector field in the adapted basis p.lOp will behave as the components of a vector 
field on the base manifold M. 

Lemma 3.6. Consider a semispray S, a nonlinear connection {ha,Vk), and the dynamical covari- 
ant derivative V associated to the pair (S, {ha, Vk))- A vector field X on T^M \ {0} is a newtonoid 
if and only if 

(3.37) x = x^'+l-y^^x^' +'-y^^x^' +... + ly-^x^ ' 



5x^ 1! 2! (5?/(2)* fc! 

where for a G {1, fc} we denote 

(3.38) 1V(")X^ = ^^"(^'') + E(^^^^(«.^"''(^^)- 

For a vector field X e j''-^X{T^M \ {0}) we have that X £ j''-^ (X|) if and only if 

Vk{X) = J{VX). 

Proof. Consider a vector field X G X{T^M \ {0}). In view of local expression p.32p and using 
notation (|3.38p . we obtain that X is a newtonoid of order k if and only if for a G {1, fc} we have 

hc,X = ^V(")X' ^ 



which is equivalent to formula p.37p . 

For a vector field X e J''^^ (X|) let Z € X| be such that X = J^'^Z. Using formula ((XT5)l it 
follows that 

J{VX) = [J oV o J''-'^)Z = [J ohk^io Csohk-^io J^-^)Z 
= (jo/ifc_io/:so(j'^-l-t;fcO J^-1))Z. 

Using J o hk-i = Wfe o J, the fact that for Z G X| we have (J o Cs o j''^'^)Z = (1 - k)J^^'^Z, and 
the second formula in p.28p we obtain J{S/X) = Vk{X). 

Suppose X = J^-^Z for some Z G X{T^M \ {0}). Starting with Vk[X) = J {SIX) then the 
arguments used to prove the first implication show that VkJ^CsZ = 0. Then 

J^-^-K^Z = J^-^{ld+JCs)Z = J^-^Z + J^CsZ = X, 

and X G J'^^i Im7r| = J^-^X%. □ 
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For a — 1, formula (|3.38p reduces to formula (|3.20p . which represents the action of the dynamical 
covariant derivative V on the horizontal components of a vector field. Using formula (j3.37p it follows 
that V(")X% which are given by formula p.38p , represent the components of a newtonoid vector 
field of order k with respect to the basis (|3.10p . Therefore V'^'^^X', behave as the components of a 
vector field on the base manifold, for arbitrary components X' that behave as the components of 
a vector field on the base manifold. This observation will allow us to extend the action of V*-"-* to 
arbitrary vector fields on T'^AI \ {0} as follows. Consider a vector field X on T'^M \ {0} expressed 
in the basis p.lOp as follows 

13=0 ^ 

Hence, for each G {0,1, ...,fc}, X^^^^^ behave as the components of a vector field on the base 
manifold. Accordingly, for each a G {l,...,fc}, the components V'^-'X*^'^)' will have the same 
property. Therefore, we can define the a-th dynamical covariant derivative of X by 

(3.40) V(")X = J2 ^^"^^^''^'tW' " ^ 

Denote V(°) = Id. For each a £ {1, k}, we set V^")/ = S'"(/), for / £ C°^{T''M \ {0}). Then, 
V^"^ satisfies a Leibniz-type rule of order a: 



[3=0 



Using the Leibniz-type rule of order a, and the requirement that V*-"-* commutes with tensor 
contraction, we can extend the action of V*-"-* to arbitrary tensor fields. For example, consider w 
a 1-form on T*^M \ {0}. Then, for an arbitrary vector field X on T'^M \ {0}, we have 

(3.41) (V(")c.) {X) = 5"(.(X)) - ± (") (v("-«.) (V(«X) . 

(3=1 ^' 

We emphasize that, for the moment, V^"^ may not be the a-th iteration of the dynamical covariant 
derivative V, unless a specific nonlinear connection will be chosen. In Theorem l3.10l we characterize 
when V(") = V" for all a > 0. 

Proposition 3.7. Consider the dynamical covariant derivative V associated to a pair {S, {ha, Vk))- 
A vector field X G X{T^M \ {0}) is a dynamical symmetry if and only if it is a newtonoid and 

(3.42) (j'^V('^)x) +$(X) = 0. 

Proof. A vector field X G X{T^M\{Q}) is a dynamical symmetry if and only if ho\S, X] = 0, for all 
a G {0, fc — 1}, and Vk[S, X] = 0. The first set of conditions ha[S,X] = for all a G {0, fc — 1} 
(or the first k conditions) are equivalent with [iS, X] G Ker J, which is equivalent with X G X|. 
For any X G X{T^M \ {0}) we have 

V{vuX) + ^{X)=v,,[S,X]. 

For any for X G X| (with h^X = X^S/Sx^), formulae (P37|) and (piiO)) imply that 

Using the above observations the result follows. □ 
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From formula (j3.19p we see that the dynamical covariant derivative V has a different action on 
distributions Hq, . . . , Hk-i than on distribution Vk- The next proposition characterize when these 
actions coincide. For example, one sufficient condition is that V J = 0. When fc = 1 we know that 
the canonical nonlinear connection for a semispray is uniquely determined by V J = 0. Proposition 
13.81 shows that when fc > 1 there is more freedom, and VJ = determines only part of a nonlinear 
connection. To uniquely determine a non-linear connection from a semispray one needs stronger 
conditions. Such conditions will be given in Theorem 13. 101 

Proposition 3.8. Consider a semispray of order k, S, a nonlinear connection {ha,Vk), and the 
dynamical covariant derivative V associated to the pair (S, (ha, Ufc)). Then, the following conditions 
are equivalent: 

i) V restricts to a map V : j''^^{Xg) — >■ j''^^{Xg) that satisfies the Leibniz rule with respect 
to the * product; 

ii) VJ" = 0, for some a G {1, fc}; 

iii) (A; + l)vk o J''-^ = J''-^ + CsJ o J^'^ ; 

iv) Nl^^=dG'/dy^^^^- 

v) V = V^^\ where V^^^ is defined in formula (|3.40p . 

Proof. We first prove that V = if and only if V J" = 0, for some a ^ {1, ...,k} . In view of first 
formula (|3.24p it follows that if V J" =0, for some a £ then V j'^ =0. For the converse 

we will show that V J*"' ~ implies V J = 0, which will also imply V J" = 0, for all a S {1, k}. 
From first formula p.24p . we have that VJ*"' = VJo J^^^ and therefore the (1, l)-type tensor field 
VJ vanishes on vector fields in Im j'^"^ = Vk~i- From local formulae p.25p it follows that VJ 
vanishes on vector fields in ®^Zq Ha, and hence V J = 0. 

Using the above considerations, formulae p.24p and (|3.25p . and the identity CsJ^ = CsJ o 
J^^^ — [k — \)J^~'^ we have that conditions ii), iii), and iv) are equivalent. 

We prove implication ii) i). Consider X S J^~^ (-^s)- According to Lemma [3.61 we have 
that VkX = JVX. If we apply V to both sides of this equality and use Vwfc = and VJ = 
we obtain Vfc(VX) = JV(VX), which by Lemma l3Jl implies that \7X G J''^^ Hence, V 

preserves the set J^~^{Xg). 

For / G C°°{T''M\{0}) and X G J''~^ (X|), formula (|336)) reduces to f * X = fX + S{f)JX. 
Using the fact that V satisfies the Leibniz rule for the usual • product it follows that 

(3.43) V(/ * X) - V(/) *X-f*VX = V(/) (V J) X. 

Since V J = it follows that V satisfies the Leibniz rule for the * product. 

We prove now the implication i) =^ ii). By formula (|3.25p . which gives the local expression 
for VJ", we only need to show that VJ = on vector fields in Im J*^"^. A set of generators 
for vector fields in Imj''-^ is given by j''^'^ (X|) U (X"''(r''M \ {0})). Since VJ vanishes on 
(X"*" (r'^M \ {0})), it remains to show that VJ vanishes on j'^'^ {^s)- Since V satisfies the 
Leibniz rule for the * product, using formula (|3.43p it follows that S{f) (V J) X = 0, for an 
arbitrary function / G C°°{T^M \ {0}) and an arbitrary vector field X G J^^^ (-^s)- Therefore, 
VJ = on the set j''-^ (X|). 

Using formulae ((XTO)) . (ESS]), and (piiO)) we have that V'^^ = V if and only if VJ = 0. Hence 
conditions ii) and v) are equivalent. □ 

Up to this point, the structure of Section [3] follows the structure of Section [2l and each result for 
a system of SODE has a counterpart for a system of HODE. Due to the complications imposed by 
the geometry of a system of HODE we will have results that are specific only to the higher order 
case. 
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Although Proposition 13.81 fixes only a part of the nonlinear connection, we can obtain now an 
explicit covariant form for the equations of variation (|3.3ip . given by the following Jacobi-type 
equations p.44p . 

Proposition 3.9. Consider V the dynamical covariant derivative associated to a pair (S, {h^, Vk)), 
and suppose that VJ = 0. Then, a vector field X on X[T^M \ {0}) is a dynamical symmetry if 
and only if X is a newtonoid and 

(3.44) Iv (V(^)X') + J2 ^^(.b V^")^^" = 0, 

0=0 

where X^ are the horizontal components of X such that hoX ~ X^5/Sx^. 
Proof. If X is a newtonoid, then formulae p. 141) and (|3.37p imply that 

and the result follows by Proposition 13.71 □ 

Note that while using formula p.42[) in Proposition 13. 7i to characterize a newtonoid X on 
T'^M \ {0}, we essentially need the assumption V J = 0. This implies that the two actions of V in 
formulae p.l9p coincide. Therefore, the action of V on the components of the vertical vector field 
VkX is given by the second formula (|3.20p . 

Next theorem will completely determine the nonlinear connection and this will allow us to 
provide simpler expressions for the k curvature components, of the Jacobi endomorphism 

Theorem 3.10. Consider a semispray of order k, S, a nonlinear connection {ha,Vk), and the 
dynamical covariant derivative V associated to the pair (S, (ha, Vk))- Then, the following conditions 
are equivalent: 

i) CsJ + ld-{k + l)vk = ij'^; 

ii) V(") = V", for all a 6 {1, 2, k}; 

in) Mi^^^ = 9GV5y(^-)' and aAq^^^ = S (mI^^,^,) + Mf^.^j^-Mj^^^, for all a e {2, k}; 
iv) Ni^y^ = d&/dy<^^)^ and aiV^^^,^. = S - Nl^_^^^Nl^^, for all a e {2, k} . 

Proof. Using formulae p.l4p and p.26p we obtain that equality i) holds true if and only if V J = 
and for each aG{l,...,fc — l}we have 

In view of formulae (|3.17p , for a in { 1 , . . . , A; — 1 } , above equations (|3.45p are equivalent to 

k — a 

(3.46) {k + 1- a)Nl,^,_^^^ = S{Nl,_^^^) - N(,^^Nl,_^)^ - ^ /f^),iVi,+i_„_^) . 

When a € {l,...,fc — 1}, equality (|3.46p is equivalent with = 0. Hence condition i) 

holds if and only if V J = and P^^^-^j — for all a € {2, fc}. Induction and equation ()3.15p show 
that the latter conditions are equivalent with the second formulae in iv). 
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Let us first note that condition ii) is equivalent with VJ = and V'"-* ~ V o V^" ^' for all 
a G {2, . . . , fc}. For a e {1, ... , k}, and /3 £ {0, k}, equations ((X^ and imply that 

s s 

and when a € {2, . . . , fc} we have 

It follows that ii) and iii) are equivalent. 

Using formulae p.l2p and (|3.4ip . and induction for a G {1, ■■■,k}, it follows that V'"-* has the 
following action on the dual basis p.lip 

V(")jy(«^ = -a!iV('„)/2/(«^ /3 G {0, 1, fc}. 

From this formula, using the action of the dynamical covariant derivative V on 1-forms, we obtain, 
for each G {0, 1, fc}, 

(V o V("-1)) = (« - 1)! [S (7V('_,)^.) - A^(^„„i)piVf^) J Si^^^. 

In view of the above two formulae it follows that condition ii) for a G {2,...,fc} and the second 
formulae in condition iv) are equivalent. □ 

The (1, l)-type tensor CsJ — Id+(fc + 1)?;^ in the first item of Theorem 13.101 can be used to 
characterize some other nonlinear connections. For example if we require it to vanish, we obtain a 
nonlinear connection studied in [U [6l [181 HD] . 

For a semispray S of order fc, we say that the nonlinear connection characterized by any of the 
four conditions in Theorem l3.10l is called the canonical nonlinear connection of S. For the remaining 
of this section we will consider this canonical nonlinear connection induced by a semispray S. The 
dual coefficients MJ^^s^^ of this connection were considered first by Miron and Atanasiu in [3^1 . Next 
theorem characterizes the symmetries of a semispray S using the dynamical covariant derivative 
V associated to a semispray S and its canonical nonlinear connection. 

Theorem 3.11. Let {ha,Vk) be the canonical nonlinear connection induced by a semispray S, and 
let V be the corresponding dynamical covariant derivative. 

i) A vector field X on T'^M \ {0} is a dynamical symmetry if and only if X is a newtonoid 
and 

(3.47) -?-V'=+ + 4>(X) = 0. 

fc! 

ii) A vector field X G X(A/) is a Lie symmetry if and only if 

(3.48) _Lv''^+^X'='° + ^iX^'^) = 0. 

fc! 

iii) Locally, a newtonoid vector field X on T^M \ {0} is a dynamical symmetry of S if and 
only if its horizontal components satisfy the Jacobi equation 

k-i 

(3.49) ry V^^'+^X^ + J2 —^Ui^"^' = 0- 



a=0 
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iv) The components of the Jacobi endomorphism are given by 



(3.50) Rl)^ = ^k + l)i—--Nl,^,_^-^^], ae{l,...,k-l} 



Proof. Part i) follows by combining Proposition 13.71 and Theorem 13. 101 Part ii) follows from Part 
i) using the fact that J^X^'^ = X^'^ . Part iii) represents the local expression of Part i). In part 
iv) , the first formula in p.50p was proven in the proof of Theorem 13.101 The second formula in 
p.50p follows by setting q = and P^^-^. = 0, . . . , = in formula p.l7p . □ 

3.3. First order variation for systems of HODE. Let c : / — !> M be a geodesic of a semispray 
S of order k. A geodesic variation of c is a smooth map V : I x {—e,e) — > M , V = V{t,s) such 
that 

i) V{t,0) = c{t), for all t in /, 

ii) V{t,s) is a geodesic for all s in (— e,e). 

Definition 3.12. A vector field ^ : / — ;> TM, along a geodesic c = tt o ^ of a semispray S, is called 
a Jacobi field of S if it is the variation vector field ^(t) = dsV{t, s)|j,^q of a geodesic variation V. 

Next we provide a sufficient condition for a vector field ^: / — > TM to be a Jacobi field. This 
sufficient condition generalizes the traditional Jacobi equation to a semispray of order k. 

For a vector field ^ : / — > TM, ^(i) = (a;'(t), ^'(i)) , along a regular curve c = tt o ^, and 
a e {0, k}, we define the lifted vector fields : / T{T^M \ {0}), along the fc-th jet lift j'=c 
of c, as follows 

Since ^ is a vector field along c, for each t G I, we can find a vector field X G X{M) such that 
X oc = ^ near t. Then ^'^■^ = X'^'" oj^c, for all a £ {0, fc}. We assume now that c is a geodesic, 
which means that j'^'c is an integral curve of a semispray S. For a > we define 

(3.51) ^ (V"X'^''") o j'^'c. 

If locally 5 = ix\C) and X = X'd/dx\ then V^'^^o = V S,' d / dy'^^^', where 

By induction on a we find that 

V"f = (V"X*) o /c, a > 0, 

where V"X* are iterations of formula p.20p . that correspond to formulae p.38p for the nonlinear 
connection determined by Theorem 13. 101 Hence, using V"X*^'° = V"X*9/9t/'*''-'', 



It follows that definition p.5ip does not depend on the choice of X, and components V"^' transform 
as a tensor on M. 

Proposition 3.13. Let ^ be a vector field ^ : I ^ TM along a geodesic of a semispray S of order 
k. If ^ is a Jacobi field of S then it satisfies the following Jacobi equation 

(3.53) l^yk+i^k,o ^^(^^k,k^^ 

fc! 
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If locally ^(t) = (a:;*(t), ^*(t)), equation (|3.53p is equivalent to 

a— ^ 



Proof. By induction on a we have that 

(a)i/:k,k 



1 



^j^tai^^/o^/c ^ — v"r, ae{0,...,fc}. 



Hence formulae p.53p and p.54p are equivalent. To prove that formula (|3.53|) holds, let c: I ^ M 
be the geodesic c = tt o ^ of 5, let t g /, and let X be a vector field X G X(M) such that X oc = ^ 
on some interval /q 3 t. Since ^ is a Jacobi field there is a geodesic variation V{t, s) such that 
V{t,Q) = c{t), dV'/ds{t,0) = Cit), and 



Differentiating the last equation with respect to s and setting s = gives 

Since c is a geodesic, it follows that S'"(/) o j''(c) = d" {f o j'' c) / dt" for all a > and / G 
C°°(r'=Af \ {0}). Thus 

(S''=+i(X') + (A: + iy.X''-^''iG')) o j'^c = 0. 
By Lemma [33] and inclusion (|335l) it follows that J[S,X'''''] = 0. Then 

Vk[S,X''-'''] = ^S''+\X')+X''^''{G'), 

and [S,X'^''^] o j'^c = 0. Repeating the argument in Proposition 13.71 shows that 

Equation (|3.53p follows. □ 

A proof for the converse of Proposition 13.131 will require to extend all the techniques developed 
in [11]. Another aspect that will have to be addressed in the future is the role of the curvature 
components R\^^)j for the geodesic behavior. 

4. Applications 

In this section we motivate the applicability of the geometric theory developed in this paper 
using examples from various fields such as: the equivalence problem, the inverse problem of the 
calculus of variations, and biharmonicity. 

First, we consider a second order Lagrangian, L2, derived from a Riemannian structure. The 
corresponding Euler-Lagrange equations form a system of fourth order differential equations, whose 
solutions are biharmonic curves [13]. For this system of differential equations, the components of the 
Jacobi endomorphism are functions of the curvature of the Riemannian metric and its dynamical 
covariant derivatives. This will motivate the use of curvature components for the components of 
the Jacobi endomorphism. 

We also show that the curvature components of the Jacobi endomorphism are useful to express 
geometric invariants that were associated previously to third or fourth order ordinary differential 
equations. See [20l EB ESI [31] . 
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4.1. Prolongation of a Riemannian structure. In this section we start witli a Riemannian 
space {M,g) and construct a 2-nd order Lagrangian L2 G C°°{T^M) [3]. Tire Euler-Lagrange 
equations for L2 determine a system of fourtlr order ordinary differential equations and hence a 
semispray of order 3. Solutions of this system are biharmonic curves for the Riemannian space |13) . 
The 3 components -R(Q,)j, G {0, 1,2}, of the Jacobi endomorphism can be expressed in terms of 
the curvature components Rji^i and their first and second order dynamical covariant derivatives. 

Consider g = gij{x)dx^®dx^ a semi- Riemannian metric on M and denote by 7jfc(a;) its Christoffel 
symbols. The geodesic spray of the Riemannian metric g has the following components of the Jacobi 
endomorphism i?^(a;, y^^-*) = R\^i{x)y'^^^^ y'^'^''^ . We denote 

(4.1) z(2)^ = y(2)^ + l7,\.(x)y(^'^'2/^^^^ 

It follows that z*^^)* behave as the components of a vector field on M. These components were 
interpreted as the covariant form of acceleration in |9l (6.5)] and as half of the components of the 
tension field in [13]. Therefore, the function L2 : T^M R, given by 

(4.2) L2(x,j/(i),y(2)) = i.9.,z(2)^z(2)^ 

is a second order Lagrangian, (U [13] . The variational problem for L2 leads to the following Euler- 
Lagrange equations 

(43) aL,_ d Id^/^X^ 

which form a system of fourth order ordinary differential equations. Consider S G X{T^M) 
(4.4) S = yi^^^l-+2y^^^^' +3yi^^^J---4G^ ' 



the corresponding semispray of order 3. This semispray is uniquely determined by the following 
equation 

dL2 ^ ( dL2 \ , I r,2 f ^-^2 



Indeed, the functions G", of the semispray S, uniquely determined by equation (j4.5p . are given by 
M «^ , f dL2\^ 1,2 ( dL2 

where dx is the Tulczyjew operator on T^M . Although, above expression for the coefficients G* of 
the semispray S is not so easy to handle with, we will show that one can obtain the corresponding 
nonlinear connection, which is determined by any of the four equivalent conditions in Theorem 
WM In formula we apply d/dy'-^^'', use the fact that the Lagrangian L2 and its partial 

derivatives are functions on T'^AI we obtain 

^ A^(\),(x,2/(^') = Mi,)^.(x,,(i)) = 7^,(x)yW'=. 

For the other two coefficients N^^^^{x, y'^-*, y*-^-*), TV^g^^ (x, y'^^\y'-^\y'^^^), of the nonlinear connection 
we use formulae iv) from Theorem 13. 101 

Consider V the dynamical covariant derivative associated to the semispray S and the corre- 
sponding nonlinear connection. It follows that 

z(2)^ = ivy(i)\ 
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hence the name of covariant form of acceleration for z'^'*. Equation (14. 5p is then equivalent to 



SL2 f SL2 \ 1 „2 ^^'2 



Using the following formulae 

we can rewrite equation (j4.7|) as follows 

(4.8) V2z(2)» + i?jz(2b' = 0. 

Note that this equation is equivalent to 

Therefore the Euler-Lagrange equations (|4.3I) can be written as follows 

which is the system of differential equations for biharmonic curves [131 (I)]- Since a geodesic of 
the Riemannian metric g is a solution of the system V{dx'^/dt) = 0, it follows that any geodesic 
is a biharmonic curve. Conversely, if the sectional curvature is non-positive then any biharmonic 
closed curve is a geodesic, see [55]. However, if we change the parameter of an arbitrary geodesic 
by a third order polynomial, the new curve is no longer a geodesic but a biharmonic curve, see [5]. 

Using any of the equivalent form (|4.8p , (|4.6p , or (|4.7p , of equations (|4.5p , we can express the 
components ^ ^ ^J" Jacobi endomorphism in terms of the curvature tensor R^jj^i 

of the Riemannian metric. Two of them are 

The other component ^(g) ^ , has a more complicated formula, but it is a function of the curvature 
tensor -R^^.;, its first, and second order dynamical derivative. 

4.2. The Wuenschmann invariant. A geometric invariant that relates third order ordinary 
differential equations with certain classes of conformal Lorentz metrics on three dimensional mani- 
folds, was proposed by K. Wuenschmann in his Ph.D thesis of 1905. The role of the Wuenschmann 
invariant for the geometry of a third order differential equation under contact transformations was 
discussed in [35]. See also for a generalization of the problem to fourth order ordinary differ- 
ential equations, where the equivalence problem leads to three invariants. A geometric method for 
obtaining the Wuenschmann invariant for a system of third order ordinary differential equations 
was proposed in |20) . In this section we show that the Wuenschmann invariant can be obtained 
from the components of the Jacobi endomorphism. 

On a 1-dimensional manifold M, consider a third-order differential equation 

and the corresponding semispray S of order 2. There exists an associated conformal Lorentzian 
structure on the 3-dimensional solution space of this equation if and only if an invariant vanishes, 
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This invariant is called the Wuenschmann invariant and it is given by [381 eq. (4)] 



^n2f 9G \ ^ dG ^ f dG \ ^ f dG \ ^ f dG \^ ^ ^ dG dG ^dG 



Consider the two curvature components given by formula p.50p of the Jacobi endomorphism: 



dG ^ dG dG ^c.2fdG\^fdG 



^ dx ^ dy^^) dy(.^) 2^ ydy(.^) J ^ {dy( 
dG 3 f dG V dG 
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dyW 2 Vay(2) J 2 ydyi^) 

It follows that the Wuenschmann invariant and the components of the Jacobi endomorphism are 
related as follows 

(4.9) = Vi?(i) -2i?(o). 



It is important to note that same formula (|4.9p was obtained in [2^ for two other different nonlinear 
connections. 

4.3. The inverse problem of Lagrangian mechanics for scalar fourth-order ordinary 
differential equations. In [22], Fels has shown that a scalar fourth-order ordinary differential 
equation admits a variational multiplier if and only if two invariants vanish. These two geometric 
invariants were associated with a fourth-order equation using Cartan's equivalence method. One 
of these invariants appears also in the list of invariants proposed by Dridi and Neut in [5T] for 
studying the equivalence problem for fourth order differential equations under fiber preserving 
diffeomorphisms. In this section we show that this invariant can be expressed in terms of the 
components of the Jacobi endomorphism using a very similar formula as in the previous section. 
On a 1-dimensional manifold M, consider a fourth-order differential equation 

d'^x ( I dx I d^x I d^x\ 

— +4.G yx,- — — — j -0, 

and the corresponding semispray S of order 3. Consider the following Wuenschmann- type invariant 
2^^ f dG \ dG / dG \ / \ ^ f dG V , „ dG dG „ dG 



3^ [dy(^)) '^dy(3-)^[dy(^->)^'^^{dy(^)) 3 [dyi^) ) ^ ^ dy(^) dyi^) ^^yd) 

The invariant W4 is related to the invariant Ii considered by Fels in [22] as follows W4 = — 3/i. 
Also (8/3)11^4 is the numerator of /g in [21] . 

Consider two curvature components and i?(2) given by formula p.50p of the Jacobi endo- 
morphism: 

R _ . dG dG 2 f dG\ 4/ dG \' 

^(1) *5y(i) 9y(2)9y(3) 3 {dy^^^ ) 3 \dy(3) J ' 

dG dG V dG 



R(2) = - 2 I 1 - 25 



5y(2) \dy^^) J \dy(^'> 

It follows that the invariant W4 and the components of the Jacobi endomorphism are related as 
follows 

(4.10) W4 = Vi?(2) -2i?(i). 

We note the similarities for expressing the two invariants W3 and W4 in formulae (|4.9p and (|4.10p . 
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As an example consider the classical spinning particle |31) . The motion of a particle rotating 
about a translating center is governed by the following system of fourth order differential equations 

(4.11) ^ + = 

where w is a real, non-zero constant. For this system of fourth order we have 
The only non-vanishing component of the Jacobi endomorphisni is 

Therefore, the Wuenschmann invariant is W4 — Vi?(2) ~ = 0. Note that the system (|4.1ip is 

separable and hence we can view each of the three equations of the system separately. For each of 
them, the Wuenschmann invariant vanishes and therefore each equation is a Lagrangian equations. 
A second order Lagrangian for the system (|4.1H) is the following 

For the system (|4.1ip . the Jacobi equations are 

V^f + w^v^f = 0, 

which can be integrated. Therefore we can obtain all geodesic variations of the system (|4.1ip . 
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